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Abstract 

We discuss a general, exact (in that matrix elements are preserved) mapping 
of fermion pairs to bosons, and find a simple factorization of the boson rep- 
resentation of fermion operators. This leads to boson Hamiltonians that are 
Hermitian and finite, with no more than two-body operators if the fermion 
Hamiltonian has at most two-body operators, and one-body boson transition 

operators if the fermion transition operator is one-body. 
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Pairwise correlations are often important in describing the physics of many-fermion sys- 
tems. The classic paradigm is the BCS theory of superconductivity where the wavefunc- 
tion is dominated by Cooper pairs which have electrons coupled up to zero linear momentum 
and spin; these boson-like pairs condense into a coherent wavefunction. Another example is 
the phenomenological Interacting Boson Model (IBM) for nuclei 0] , where many states and 
transition amplitudes are successfully described using only s- and d- (angular momentum 
J = 0, 2) bosons, which represent coherent nucleon pairs. In both cases the large number of 
fermion degrees of freedom are well modelled by only a few boson degrees of freedom. On 
the other hand, however, despite some forays by Otsuka et al. 0, a rigorous microscopic 
basis for such phenomenological models is lacking. 

The basic problem is to represent the underlying fermion dynamics and statistics with 
a boson image amenable to approximation and numerical calculation. Considerable effort 
has gone into mapping fermion pairs into bosons 0,^. However these mappings typically 
suffer from a variety of defects. Most, such as the Belyaev-Zelevinskii ^ and Marumori ^ 
mappings give rise to boson Hamiltonians with infinite expansions, that is, A^-body terms 
with oo. Convergence is slow even when "collective" fermion pairs are used. Finite 

but non-Hermitian boson Hamiltonians have also been derived [§,0. 

In this Letter, using an alternate approach, we show that the infinite expansion boson 
Hamiltonian obtained from the exact mapping of fermion matrix elements to boson matrix 
elements can be factorized into a finite, Hermitian boson Hamiltonian times a norm operator, 
and it is the norm operator which has an infinite boson expansion. 

Consider a fermion space with 2Q single-particle states, and a fermion Hamiltonian H. 
The general problem is to solve the fermion eigenvalue equation 

H\^p)=E,\^f,), (1) 

find transition amplitudes between eigenstates, and so on. To do this we require a many- 
body basis. Often the basis set for many-fermion wavefunctions are Slater determinants, 
antisymmeterized products of single-fermion wavefunctions which we can write as products 
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of the Fock creation operators Oj, j = 1, ■ ■ -,2^1 on the vacuum aj^ ■ ■ ■ aj^^ |0) for n fermions. 
These states span the antisymmetric irreducible representation of the unitary group in 2Q 
dimensions, SU(2fi). But for an even number of fermions one can just as well construct 
states from N = n/2 fermion pair creation operators, 
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M = n 4™ io) ' (2) 

m=l 

with 

(3) 

We always choose the Q{2Q — 1) matrices to be antisymmetric to preserve the underlying 
fermion statistics, and we choose for the normalization the trace tr (^Aq,A^^ = 6ai3- For this 
Letter we represent generic one- and two-body operators by 

f^Y.T^Jala„ V^^(5|V|A)iii,. (4) 

ij SX 

We begin with the straightforward mapping to boson states 

N 

m - \M = n io) ' (5) 

m=l 

where the are boson creation operators. In conjunction with this simple mapping of states 
we construct boson operators that follow the philosophy of Marumori and preserve matrix 
elements, for example introducing boson operators Tb, Vb, and most importantly Tig such 
that 7b 10/3) = {ipa\f {(f)a\ Vb 10/?) = M V , and (0o| Hb 10/?) = {ipa\ H \lpp). 

In addition, because wavefunctions of the form @ do not form an ortho normal set, we 
construct the norm operator A/b with the property {4>a\-f^B \4>i3) = ("^a li'is)- While the 
mapping (^) is the usual Marumori mapping, this boson representation of fermion operators 
is different, stemming from the fact that Marumori does not have an explicit norm operator. 
With this mapping the fermion eigenvalue equation (||) becomes a generalized (because of 
the norm) boson eigvenvalue equation 

nB\%) = Epj^B\%)- (6) 



Because we have defined our boson operators so as to preserve matrix elements, the original 
energy spectrum of (|I|) is found in (^. However, because the boson space is much larger than 
the original fermion space, also has additional spurious states that do not correspond to 
physical fermion states. These by construction have zero eigenvalue. 



When one constructs the norm operator one finds it can be conveniently ex- 
pressed in terms of the kth order Casimir operators of the unitary group SU(2i7), Ck = 

2Ht {j:^,biKA^Ai)\ 



:expff;^t!c,): 



\k=2 



k 



(7) 



where the colons ':' refer to normal-ordering of the boson operators. Expanding in 
normal order one obtains the form 



= 1 + E E ({^} ; M) n n K ■ 



(8) 



i=2 o-T 171=1 n=l 

The £-body boson terms embody the fact that fermion pair creation and annihilation oper- 
ators do not have the same commutation relations as do boson operators, and act to enforce 
the Pauli principle. 

The operators 7^, are also complicated many-body operators similar in form to (^, 
though upon examination one can write them compactly using the norm operator. For 
example. 



7b=:E(-1 



2nr 



k=l 



(TAtA.) biK{Y.^aAlbib 

\al3 



(9) 



where (T)^^ = Tjj, and similarly for Vb [CT- general these boson operators do not have 



good convergence properties, so that truncation of (|^), (^) as written is problematic. 

A key result of this Letter, as suggested by the explicit form of (|^), is that these operators 
factor in a simple way: Tb = N'bTb = TbN'b and Vb = J^bVb = Vb-^b, where the factored 
operators Tb, Vb, which we term the boson images of T, V. have simple forms. For example, 
a one-body fermion operator has a one-body boson image 
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TB = 2j2tT (TAtA^jfotfe^. (10) 
To prove this factorization, one puts TbJ^b into normal order: 

2 tr (TAt A,) {blAfBbr - hi [Nb, K] } • (H) 

err 

The first term is the k = 1 term from (j^) and the second term, the commutation of the 



VB = Y.{^\y\>^) 
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norm operator with 6,-, gives rise to the k> 1 terms |jTO[- Similarly, 

+ 2 E E (a^aJA,, At A;,At,) hlhlKK' (12) 

cfcf' tt' 

and in general one can find a image Hamiltonian 

HB=fB + VB. (13) 

Most of the complexity resides in the norm operator and norm matrix and the slow conver- 
gence of many-body terms that arise in other boson mappings most likely do so because the 
norm is not completely factored out. 

While we can demonstrate this factorization using the explicit forms for 7^, Mb-, etc., 
|T0| , one can arrive at this result directly. For example, the action of a one-body operator T 
on a fermion state |'?/'/3) is to linearly replace one fermion creation operator with another, or 
in the fermion-pair representation (0) , one fermion-pair creation operator with another, and 
so schematically T linear combination of This replacement is easily represented 

by bosons and it is straightforward to show that in fact 

f 1^^) = i^a 1^,) (0.1 Tb |0;3) = E (^c^ITb 10,) m (14) 
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with Tb defined as above and without any explicit reference to the norm operator, and 
similarly for a two-body fermion operator V. 

Thus any boson representation of a Hamiltonian factorizes: "Hb = MbHb- Since the 
norm operator is a function of the SU(2fi) Casimir operators it commutes with the boson 
images of fermion operators and one can simultaneously diagonalize both Hb and Mb- 
Then Eqn. (^ becomes 



Hb\%) = Ep\%). (15) 

where Ep = Ep for the physical states, but Ep for the spurious states is no longer necessarily 
zero. 

The boson Hamiltonian Hb is by construction Hermitian and, if one starts with at most 
only two-body interactions between fermions, has at most two-body boson interactions. All 
physical eigenstates of the original fermion Hamiltonian will have counterparts in (p^f). It 
should be clear that transition amplitudes between physical eigenstates will be preserved. 
Spurious states will also exist but, since the norm operator A/^ commutes with the boson 
image Hamiltonian Hb, the physical eigenstates and the spurious states will not admix. 
The spurious states can be identified because the norm operator As annihilates spurious 
states and has eigenvalue {2N — 1)!! on physical states. To aid in the convergence to the 
physical states in numerical computations we can put them low in the spectrum by adding 
the normal-ordered second order Casimir operator minus its eigenvalue for the physical 
states, 

M = g[.C2:+4:N{N -1)] (16) 

where g is an arbitrary strength large enough to lower the physical states below the spurious 
states. Using the known eigenvalues of the Casimir operators for the irreducible represen- 
tations of SU(2fi) [|TT|, we find that M has zero eigenvalue for the physical states (those 
belonging to the totally antisymmetric irrep) and a positive definite spectrum for the spu- 
rious states. 

We have mapped the fermion Hamiltonian into a finite, Hermitian, boson Hamiltonian. 
While diagonalizing this boson Hamiltonian yields a spectrum that contains the eigenvalues 
of the original fermion Hamiltonian, such a diagonalization would be numerically more 
taxing than the original problem because of the additional spurious boson states. Instead 
one wants to restrict the problem to a few boson degrees of freedom. There are three steps 
in this process. 
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The first step is identifying tlie dominant degrees of freedom, either boson or fermion- 
pair. This may be done by assumption (prejudice), such as choosing J = 0, 2 fermion pairs 
or s- and rf-bosons to match the IBM, or by some variational method using a fermion-pair- 



condensate 121 or Hartree-Bose. 



From a restricted set of fermion pairs, one can construct wavefunctions of the form (2). 
Diagonalizing the fermion Hamiltonian in this subspace yields approximations to the true 
eigenstates; how good an approximation depends on the first step. The second step then is 
to find appropriate boson operators for the restricted space, that is operators that preserve 
the fermion matrix elements and keep physical and spurious states from mixing, using only 
those bosons corresponding to the allowed fermion pairs. Construction of these operators 
cannot be carried out naively. The boson representations of the norm operator and of the 



and 



n 



B 



respectively, are 



Hamiltonian in this restricted space, which we denote by Mb 
easily found: one takes the boson representations in the full space (for the norm operator 
Eqn. (^, (^) and keep only those terms consisting of allowed bosons. The coefficients are 
unchanged and matrix elements of states in the truncated fermion space are still preserved. 



Factorizing the Hamiltonian is not so simple, however, as 
make a factorization 



H 



B 



T 



H 



B 



Me 



Ml 



B 



One can 



(17) 



but Ht 7^ Hb ^. Instead one starts from a general ansatz which we simply sketch as 



Hj, = ei¥b + e2b^¥bb + e^b^b^b^bbb + ... 



(18) 



and find the coefficients Oi (which depend on the truncation) by normal-ordering the right- 
hand side of (0) and equating term by term with the left-hand side. The result is a boson 
image of the Hamiltonian whose spectrum contains the eigenstates not of the full fermion 
space but of the restricted fermion space. 

The final step, which is beyond the scope of this Letter to discuss in detail, would be to 
renormalize the operators from the second step into effective operators which account for the 



excluded degrees of freedom, producing eigenvalues and transition amplitudes corresponding 
to the full fermion space. If the assumptions of the first step are valid these corrections are 
small. 

As an example we consider the pairing interaction. For fermions with spin, the dimension 
2Q is even and for, each fermion state a|, there is a time- reversed state, Kalx"^ = e^al {K is 
the time-reversal operator) where = 1 or —1 and Cj = — e^. The paired two-particle state 
is then Al = (4i7)~^/^ X^i Qolal and hence (Ao)jj = (2f2)~^/^ej5jj. The pairing interaction 
P = AqAq maps to the boson image 

Pb = hlho + 2 ^ ^ tr (A.A|,A., AtAoA|,) hlhlhrK'. (19) 

The boson image Pb can be solved completely and easily because, like the original fermion 



Hamiltonian P [0, it is a linear combination of Casimir operators, 

1 



2N + C2- 2Gi 



(20) 



where G2 is the Casimir operator of the symplectic subgroup Sp(2n); it leaves the paired 



boson invariant, 



G2, bo 



0. Because the eigenstates of Pb must be simultaneously eigen- 



states of Afs, the naive guess 



(blf |0) 



(21) 



does not yield the ground state. This is not surprising since simple correspondence, whereby 
the ground state of the pairing Hamiltonian maps into (bl)^ |0), leads to non-Hermitian 
Hamiltonians |]T^. The norm operator A/^ acting on the ansatz (^) projects out the exact 
ground state; that projected wavefunction however contains other bosons. 

If our truncated space consists of only fermion wavefunctions constructed from Aq pairs, 
and so in the boson space we have only fej-bosons, then 



B 



T 



^obo + ^blblbobo. 



(22) 



This has (^) as its only eigenstate but yields the wrong ground state energy, and one finds 
explicitly that 
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V 



B 



B 



(23) 



Instead, the boson representation 



V 



B 



correctly factorizes to Pt 



Ml 



B 



T 



with 



Pt = blbo - ^blblbobo 



(24) 



This boson image, which is found starting from the ansatz (|TH]), does indeed yield the correct 
ground state energy with ( ^l]) as its eigenstate. 

The factorization of a boson representation of a Hamiltonian T-is is guaranteed to yield a 
finite boson image only in the full space; an arbitrary truncation with a general Hamiltonian 
can induce higher-body terms, as allowed for in the ansatz (|T8D. The important question of 
the convergence of the many-body terms induced by the truncation will be explored in the 
near future. 

In summary, the mapping of a fermion Hamiltonian into a boson Hamiltonian , where the 
bosons represent correlated fermion pairs, factorizes into a finite boson Hamiltonian times 
a boson norm operator with an infinite number of boson terms. Because this finite boson 
Hamiltonian commutes with the boson norm operator, the boson Hamiltonian does not mix 
physical and spurious states. 

This work was supported by the U. S. Department of Energy. 
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